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New class instabilities is identified in Hall plasmas in configurations with open magnetic field lines. It is 
shown that sheath resistivity results in a robust instability driven by the equilibrium electric field. It is conjec- 
tured that these instabilities play a crucial role in anomalous transport in Hall plasmas devices. 



1. Introduction. Plasmas in contact with material sur- 
face form a narrow region of noncompensated positive charge, 
or sheath, which occurs as a result of fast losses of neg- 
ative charge due to the high mobility of electrons. This 
basic phenomena is of fundamental importance for various 
plasma physics applications. The width of the sheath re- 
gion, typically of the order of the electron Debye length, 
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Ad = (Te/47re^rio) , is normally many orders of magni- 
tude smaller than the characteristic length scale of the bulk 
plasma region. Despite such large length scale separation 
sheath may critically affect the properties of bulk plasma, 
e.g. anomalous (turbulent) transport in magnetically confined 
plasmas HH. 

While the turbulent transport in plasmas has been studied 
for many years and there have been significant advances in 
understanding of dynamics of drift wave turbulence in mag- 
netized plasmas, the complete understanding remains elusive. 
One of the less studied field is the anomalous mobility of elec- 
trons in Hall plasmas PO. Hall plasmas, with a characteristic 
geometrical length scale between the ion and electron Lar- 
mor radii, pe ^ L ^ p,, and intermediate range of the 
frequencies between ion and electron cyclotron frequencies, 
LJci -^ u! -^ Wee, represent a special, but common and im- 
portant case of magnetically controlled plasma. Hall plasmas 
conditions are typical for many laboratory and space plasmas, 
e.g. for some regions in Earth ionosphere and magnetosphere. 

Hall plasma conditions are typical for many technological 
applications. Common feature of these devices is presence of 
stationary externally applied electric field Eq, which is per- 
pendicular to the equilibrium magnetic field Bq, and produce 
stationary Eq x Bq drift velocity. The ions, due their large 
Larmor radius, are unmagnetized and accelerated in the Eq 
direction, while the electron collisions lead to a finite current 
along the Eq. As a result, quasineutral plasma is accelerated 
along the Eq. Such plasma accelerators, typically in coax- 
ial geometry with radial magnetic field Bq, axial Eq, and 
azimuthal drift Eq x Bq , also used as technological plasma 
sources, have recently became a subject of growing interest 
due their applications in electric space propulsion, so called 
Hall plasma thrusters. 



Despite many years of successful operation in space ||4|], 
the physics of the Hall thrusters is not completely understood. 
Ongoing expansion of applications of Hall thrusters for var- 
ious space missions of growing complexities and associated 
new requirements dictate the need for better understanding 
of physical processes in these devices 151]. One of the criti- 
cal problems, is the nature of the electron current across the 
magnetic field. Electrons are nearly collisionless in typical 
applications, and coUisional values of the current are too low 
to explain the experimental data. It has long been known 
that Hall thruster plasmas are typically in a turbulent state 
exhibiting a multitude of instabilities and fluctuations in the 
wide range of frequencies [6>]. Plasma turbulence was sug- 
gested to be the reason for anomalous electron mobility, how- 
ever the exact nature of responsible instabilities remains un- 
clear At the same time, an alternative mechanism of elec- 
tron conductivity created by electron collisions with the walls, 
so called near-wall conductivity, was proposedf4j- However, 
the reflection of electrons from the ideal sheath would be re- 
versible, so the stochastisation mechanism of near-wall con- 
ductivity is unclear In the original workf4||, it was suggested 
that sheath fluctuations would create the required stochastic 
motion of electron, though the reasons for such fluctuations 
are unknown. Experimental data and numerical modeling 
suggest [Ty that a combination of both turbulent mobility and 
near-wall conductivity is required to satisfactory describe ex- 
perimental behavior 

In this paper, we propose a novel instability mechanism 
based on coupling of sound waves in bulk plasma with sheath 
fluctuations. Sheath dynamics, included as a boundary condi- 
tion for main plasma, provides the positive feedback rendering 
waves unstable. Such mechanism affects both turbulent mo- 
bility of bulk electrons and near-wall conductivity of electrons 
reflected by fluctuating sheath. Recent numerical simulations 
have also suggested the presence of sheath fluctuations [8, 9j|. 
There has been a number of previous studies which concen- 
trated on instabilities in crossed Eq x Bq fields Hall plasmas, 
which were caused by the equilibrium electron flow in com- 
bination with density and magnetic field gradients, resistive 
effects, and electron cyclotron effects. These instabilities may 



also be affected by sheath impedance effects and eb operative 
in addition to the instabiUty of basic sound waves whcih is 
studied in this paper. 

It was noted earlier that boundary conditions at the sheath 
amounts to a specific dissipation for bulk plasma, the so called 
sheath resistivitv lllSll . Sheath resistivity was shown to lead to 
some instabilities of magnetically confined plasmas for fusion 
applications 121 ll6(] but was not studied in Hall plasmas con- 
ditions. Original concept of sheath resistivity was developed 
for conducting walls. For such a wall, the electron and ion 
currents into the wall should not be balanced locally. There- 
fore, local perturbations of plasma density and potential lead 
to different ion and electron currents, thus producing a finite 
current into the sheath. For the dielectric wall, the current on 
the material surface must be zero. However, the total sheath 
impedance has a finite and large reactive( capacitive) part, so it 
is able to maintain a finite oscillating current on the quasineu- 
tral plasma side even when the current into the wall is zero. 
Recognition of this fact allows us to describe the sheath ef- 
fects as a boundary condition for bulk plasma oscillations. It 
turns out that in this formulation, the resistive part of the total 
sheath impedance leads to a robust instability in Eo x Bo Hall 
plasmas. 

Principally, the stationary Eg x Bg flow of magnetized elec- 
trons of the velocity is a powerful source of free energy in 
Hall plasma devices. The electron perturbations with phase 
velocities lower than the equilibrium flow velocity vo = 
cEo X B0/-B0 will have negative energy: the total energy 
of the perturbed system is lower than the energy of the un- 
perturbed state with equilibrium flow. The negative energy 
perturbations are prone to dissipation instabilities. Dissipation 
takes the energy away from the negative energy perturbation 
thus inducing the further growth of its amplitude. We illus- 
trate this mechanism with a simple example of the ion sound 
instability in Hall plasma with equiUbrium electron flow in 
Section II. The growth of the negative energy modes induced 
by sheath dissipation is the exactly the instability mechanism 
for fluctuations in Hall plasma plasmas studied in Section V. 
This instability affects both turbulent mobility and near wall 
conductivity. 

2. Destabilized ion sound waves in Hall plasmas. Ions are 
unmagnetized and described by the equation 
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We consider the dynamics along the magnetic field and in per- 
pendicular direction separately, so that 
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(2) 



where z is the direction along the magnetic field, z is a unit 
vector along the equilibrium magnetic field, and y is the unit 
vector in the azimuthal direction. Then ion continuity equa- 
tion takes the form 
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The electrons are magnetized and, in the lowest order, the per- 
pendicular electron velocity is given by the 
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Vj^o = cEo X b/Bo = —cEox/Boy is the equilibrium elec- 
tron drift due to the axial electric field Eqx- 

The electrons have large mobility along the magnetic field 
and have Boltzmann distribution of density 
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It is valid for w < k^VTe and is found from the electron mo- 
mentum balance along the magnetic field 
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dz dz 

electrons are assumed to be isothermal Tg = const. In the 
local approximation, Boltzmann electrons (|5]) and perturbed 
ion density from (O result in the dispersion relation for ion 
sound waves 
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Note that the equilibrium electron drift Vj^o does not result in 
the Doppler frequency shift because ions are not magnetized 
here. In this derivation, we do not utilize the information 
about the electron current along the magnetic field, though 
obviously such a current is finite. It can be found from the 
electron continuity equation 

e dz' 

where wo = k • V^o- Using quasineutrality, the ion and elec- 
tron continuity equations can be cast into the form that ex- 
plicitly shows a contribution of the current along the magnetic 
field 
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where J|| = J||j + J||e. Together with ([Hi, (|5]l, and ^, this 
equation, of course, will result in (|7|i- 

The instability mechanism of the discussed modes in es- 
sential way relies on the stationary electron flow due to the 
equilibrium electric field Eo and presence of dissipation. To 
highlight the physical mechanism of the instability we con- 
sider a simple model of the ion sound waves in Hall plasma in 
presence of electron collisions. Taking into account the elec- 
tron collisions, the parallel momentum balance becomes 
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Using this expression and the electron continuity one has for 
the density perturbations 
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Together with ion equation and quasineutrality it results in the 
dispersion relation for unstable ion sound waves 
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This instability can be interpreted as the negative energy 
sound wave, uj < luq, destabilized by dissipation. For weakly 
collisional plasmas, when Ae > L^JrrieJmi, the bulk dissi- 
pation due to interquartile collisions, becomes less important 
compared to the sheath dissipation lIlSll . where Ae is the elec- 
tron mean free path, and L is the characteristic length scale of 
the plasma slab (in the direction perpendicular to the sheath). 
In Hall thruster plasmas, sheath dissipation induces unstable 
global and small scale modes, which are considered next. 

4. Boundary conditions and sheath impedance. The sheath 
effects are considered by including the perturbed current along 
the magnetic field into the sheath. For stationary state, we 
assume the standard sheath model with Bohm ion current 



JiQ = encs, (13) 



and the electron current 
1 
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envTe exp (-ei^io/Te) . (14) 



The ion sound velocity is defined as c^ — Te/rrii. In stationary 
conditions the total current into the sheath is zero: Joe + Joi = 
0. 

Potential (f>, density h, and temperature fluctuations Te at 
the sheath boundary produce the perturbations of the electron 
and ion current into the sheath 
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In this model we do not consider the processes inside the 
sheath layer and do not specify how the finite current into the 
sheath is accommodated. The actual dynamics of this current 
will depend on the surface material and conditions, e.g. for the 
metal surface this current can be admitted into the wall, while 
for the dielectric the finite current into the sheath will lead to 
perturbation and rearrangement of the spatial charge distribu- 
tion in the sheath layer and on the surface. The characteris- 
tic frequency of sheath oscillations is of the order of the ion 
plasma frequency ujpi. Therefore, as long as the frequency of 
the perturbations is lower than ujpi and the amplitude is small. 
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the sheath perturbation will retain its standard structure (fT3] ) 
and ( fl4T ). and the current at the sheath boundary can be deter- 
mined as a small perturbation of (fT3l ) and (fT4l ). This condition 
alternatively can be viewed as the condition that the character- 
istic charge admitted into the sheath as a result of fluctuations 
J/cj is small compared to the charge of the stationary sheath 



J/uj < Q. 
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Then, the sheath will only experience a small perturbation and 
the current at the sheath boundary can be calculated from (fTSb 
and ( fT4] i. The charge is the sheath later Q ~ p^d, where the 
charge density pe is estimated as pe — (f>o/ (47rd^), where the 
sheath thickness is of the order of the Debye length, d ~ Xdc- 
Assuming that J ~ Joop/Te, one obtains the condition 
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which is equivalent to ( flSl ). 

5. Global sheath induced modes. Presence of the sheath 
at the plasma boundary imposed additional constraints on 
plasma dynamics via the conditions (fTSl l and ( fT6b . As a result 
of these constraints, the ion sound dynamics is modified and 
new modes appear In general, these modes have the eigen- 
mode structure that depends both on the z— coordinate and 
the perpendicular direction y. It is instructive first to consider 
the modes which have no structure along the magnetic field. 
Such modes can be derived within simple model equations av- 
eraged along the magnetic field lines. Since these modes have 
no profile along the magnetic field, d/dz = 0, we will call 
them the global modes. 

Consider the magnetic field line tube between z = —L and 
z = L and introduce the averaged potential and density 
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The averaged quasineutrality (current closure) equation ^ 
becomes 
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One mode can be obtained from (1231 ). by using the Boltzmann 
relation averaged along the magnetic field line tube 
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Equations ( |23] ) and ( l24b produce the dispersion equation for 
the damped mode 
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where Vsh = Cs/L is the effective "colUsion frequency" char- 
acterizing the sheath resistivity. 

However, the conditions ^ and ( l24l i is not the most gen- 
eral solution of the electron momentum balance equation ^. 
Indeed, integrating (|6]l along the magnetic field lines one ob- 
tains 



n, (z, y,t) = ^ (0 (z, y,t) + C (y, <)) ■ (26) 

The integration constant C {y, t), which may depend explic- 
itly on y and t is determined by boundary conditions for per- 
turbations at the sheath. Then the perturbed electron current 
at the boundary found from (fT6b is 






(27) 



The density perturbation induced by this current can be found 
from the electron continuity equation dSI giving 






(28) 



It is important to note that in this process, the plasma potential 
is also perturbed. This perturbation can be found from 
and I 
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Then, the perturbed plasma density written in terms of the 
electrostatic potential is 
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We have changed the notation for (/) to (/) to emphasize that this 
perturbation depends only in y and t and is independent of z. 
Alternatively, these plasma density perturbations can be 
found directly from the averaged equation (O and using ( fT6] l 
for the perturbed electron current 
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which results in the same expression (l30t . 

The perturbed ion density is found from the averaged ion 
continuity equation 
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Using quasineutrality, from Eqs. (l30b and ( |33] ), we obtain the 
following dispersion relation: 
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In absence of the equilibrium drift (ujq = 0), Eq. (??) describe 
sound waves damped by sheath resistivity 

w = -ikyC^/i^sh- (35) 

For large luq > Vg]^, it has an unstable root with 

w = ± (iujoklcl/vsh) (36) 



6. Small scale modes. In this section we consider the sheath 
instability by taking into account the eigen-mode structure in 
the direction parallel to the magnetic field. The ion continuity 
and motion equations result in the following local expression 
for ion density 
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The parallel component of the electron equation of motion is 
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(38) 



which, assuming constant electron temperature gives, and we 
assume (fi = (t){y, t) 
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Using quasineutrahty and Eqs. (l37l ) and ( [39l ), we have 
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This is an eigen mode equation for ion sound oscillation. 
Since ions are not subject of E x B, there is no Doppler 
shift in the mode frequency. The eigen mode equation ( |40t 
is solved in plane geometry —L<z<L with boundary 
conditions at z = ±L. The boundary conditions are defined 
by the current into the sheath (fTsT i and ( fTSI l 
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The ion and electron current are found from electron and ion 
equations of motion and continuity equations 



'J zi 
1 Q 

-i(a;-a;o)ne Tr^ze = 

e oz 



ie^TiQ dcj) 
rriiLU dz ' 



(42) 



(43) 



Boundary conditions (|4T]i, (l42]i, and (|43ll at z = ±L and ( |40| ) 
make a full eigenmode problem. 



General solution for (|40] | has the form 

(j) = Aex-p{ikzz) + B ex-p{—ikzz) , (44) 

where kz satisfy the equation 
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Using ( |42] ) and ( l43T l and boundary condition ( |4TI ) one gets the 
dispersion relation 
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This dispersion relation can be cast for even A = B and F^G- ^- G^wth rate of the unstable modes m Eqs. () and (). 
odd A = —B modes in the form ^^^ '^^^S^'- P^"sl shows zoomed in interspacing of even and 

odd modes. 





1- 


Lu{iO — ClIq) 


= itankzL, 


(47) 


r^z^s 

UJ 


1- 


LU{UJ — LOq) 

klcl \ 


= —i cot kzL. 


(48) 



Consider the long wavelength approximation k^L -^ 1. 
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and for odd modes with A 
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Using the relation for kz , the dispersion relation for the even 
mode in the long wavelength limit becomes 

There is no instability for the even mode in the long wave- 
length limit. 

Odd mode ( |50] | gives the equation for the complex wavector 
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1 - iuj/vsh ' 
Using it in the dispersion equation (l45T l we obtain 

w^ + [ivsh^o - k\cl) Lo - ivshk\cl == 0. (51) 
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This has the unstable root instability u = {—ivshi^o) 
for k± ^ . For typical Hall thruster parameters: plasma 
density tiq = 10^^ cm^'^, magnetic field Sq = 150 G, electric 
field Eq =200 V/cm, the electron temperature T^ = 20 eV, the 
length scale L = 1 cm, and lowest wavevector fcy=1.0 cm^^, 
the unstable mode has the frequency and growth rate of the 
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order of w^ — 7 —4 MHz. It scales as ^ (kyEo) ' . There 
is also a weakly unstable root for wq ^ corresponding to 
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FIG. 2: Real part of the frequency of unstable modes in Eqs. 

and (). The right panel shows zoomed in interspacing of 

even and odd modes. 



a weak evolution of the initial equilibrium Ill5n . The growth 
rate and real frequency of the unstable mode in|5T|are shown 
in Fig. 1. 

The small scale instabilities correspond to the roots of the 
equations Wh and ( |48] |. The numerical solutions are shown 
in Fig. 2. 

7. Summary. To summarize, we have identified novel class 
of negative energy wave instability operating in Hall plasma 
with collisionless unmagnetized ions and magnetized elec- 
trons in presence of the equilibrium electric field. The insta- 
bility is triggered by the dissipation due to the normal current 
into the sheath. Such instabilities may exist in the form of the 
symmetric global mode, Eq. (fSTT i. or small scale modes de- 
scribed by equations (l47land ( |48] |. The global mode is slowly 
varying along the magnetic field, with the characteristic length 
scale kzl/L, while small scale modes are characterized by 
fast variations along the magnetic field lines. The fluctuations 



associated with such modes will affect both turbulent mobil- 
ity of bulk electrons and near-wall conductivity of electrons 
reflected from the fluctuating sheath. 
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